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Abstract. We study cosmological consequences of a kinetic gravity braiding model, which is 
proposed as an alternative to the dark energy model. The kinetic braiding model we study is 
characterized by a parameter n, which corresponds to the original galileon cosmological model 
for n = 1. We find that the background expansion of the universe of the kinetic braiding 
model is the same as the Dvali- Turner's model, which reduces to that of the standard cold 
dark matter model with a cosmological constant (ACDM model) for n equal to infinity. We 
also find that the evolution of the linear cosmological perturbation in the kinetic braiding 
model reduces to that of the ACDM model for n = oo. Then, we focus our study on 
the growth history of the linear density perturbation as well as the spherical collapse in 
the nonlinear regime of the density perturbations, which might be important in order to 
distinguish between the kinetic braiding model and the ACDM model when n is finite. The 
theoretical prediction for the large scale structure is confronted with the multipole power 
spectrum of the luminous red galaxy sample of the Sloan Digital Sky survey. We also discuss 
future prospects of constraining the kinetic braiding model using a future redshift survey like 
the WFMOS/SuMIRe PFS survey as well as the cluster redshift distribution in the South 
Pole Telescope survey. 
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1 Introduction 

Cosmological observations of type la supernovae [1, 2], the cosmic microwave background 
anisotropies [3, 4], the large scale structure of galaxies [5, 6], and clusters of galaxies [7, 8] 
symbolize the precision cosmology of the present-day. The universe is undergoing a phase 
of an accelerated expansion, which raises the most challenging problem in the era of the 
precision cosmology, called the dark energy problem. The cosmic accelerated expansion is 
ascribed to the presence of an extra component of universe with a negative pressure [9, 10]. 
The cosmological constant fits the cosmological observational data, however, the smallness 
of the value cannot be explained naturally [11, 12]. The quantum field theory predicts the 
existence of the vacuum energy, which corresponds to the cosmological constant A. However, 
it has been known that there is the fine tuning problem i.e., the fact that the observed value 
of the vacuum energy or the cosmological constant, p\ ~ 10 _47 GeV 4 , is much smaller orders 
of magnitude than the value that we expect in the quantum field theory. 
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An alternative to explain the accelerated expansion of the present universe is to modify 
the theory of gravity at long distance, such as the scalar-tensor theory which is an extension 
of general relativity by adding non-minimally coupled scalar field to gravity [13-17], f(R) 
gravity altering the Einstein-Hilbert action by a general function of Ricci scalar R [18-24], 
and the Dvali-Gabadadze-Porrati (DGP) model developed in the context of the brane world 
scenario [25, 26]. However, these modified gravity models are not necessarily successful [27- 
31]. 

Recently, as an alternative to general relativity, there has been proposed the galileon 
gravity models, which are constructed by introducing the scalar field with the self-interaction 
whose Lagrangian is invariant in the Minkowski space-time under the Galilean symmetry 
d^cft — > d^cft + by, which keeps equation of motion at the second order differential equation 
[32-44, 46-59]. The simplest term of the self-interaction of the galileon model is (V0) 2 D0, 
which appears in the 4-dimensional effect theory of the DGP model. The galileon models 
generalized to a curved spacetime do not necessarily possess the Galilean symmetry, however, 
it is possible to keep the equation of motion at the second order differential equation. The 
galileon model in the framework of the Brans-Dicke theory has been studied in [32-38]. In 
addition to these works, the authors in [39-45] investigated the cosmology of the covariant 
galileon field. These models can lead the self-accelerated expansion of the late-time universe. 
This nature is applied to construct the inflation model [61, 63] (cf. [60, 62]). The galileon 
field equation of motion is the second order differential equation with respect to the time, 
which prevents the theory from appearing a new degree of freedom, and the perturbation 
of the theory does not suffer from ghost or instability problem. A remarkable property of 
the galileon model is the Vainshtein mechanism [66]. The self-interaction term (V</>) 2 D</> 
induces the decoupling of the galileon field (ft from gravity at small scale. This allows the 
galileon theory to recover general relativity around a high density region, which ensures the 
consistency with the solar system experiments. 

The kinetic braiding model is proposed inspired from the galileon model [46], which 
introduced the extended self-interaction term G(<j), X)\D(ft minimally coupled to gravity, where 
G(cft,X) is a function of (ft and X, which we defined X = —g ll ' / V fl (ftV u (ft/2. In the present 
paper, we consider the case G((ft,X) oc X n , where n is the parameter. In the original work 
in [46], the case of n = 1 is considered. Adopting an attractor solution for the galileon 
field, we find that the background expansion for general value of n realizes the Dvali- Turner 
model [67, 68], which is the phenomenologically modified Friedmann equation, In the Dvali- 
Turner model, the background expansion reduces to that of the ACDM model for n equal 
to infinity. We also find that the linear cosmological perturbations of the kinetic braiding 
model reduces to that of the ACDM model for n = oo. Thus, the kinetic braiding model 
connects the original galileon model and the ACDM model by the parameter n at least as 
for the background expansion and the cosmological linear perturbations. This interesting 
feature is one of the reasons why we consider this model. 

In general, when a modified gravity model shares a similar expansion history to those 
of dark energy models, there is a crucial difficulty in discriminating these models only with 
measurements of the background evolution. Therefore, it is worth examining the growth 
history of the large scale structure, because the degeneracy of the background evolution 
could be broken (e.g., [69, 70]). In the present paper, we examine the evolution of density 
perturbations in the kinetic braiding model in the linear regime as well as the nonlinear 
regime by considering the spherical collapse model, which might be important in order to 
distinguish between the kinetic braiding model with a finite n and the dark energy model as 
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well as the ACDM model. 

The paper is organized as follows: In section 2, we briefly review the kinetic braiding 
model, where the background evolution and the stability condition of the galileon field per- 
turbation are summarized. In section 3, constraints from type la supernovae observations 
and the WMAP cosmic microwave background (CMB) anisotropy experiment are presented. 
In section 4, we investigate the linear evolution of the cosmological density perturbations. 
We show that the linear cosmological perturbations of the kinetic braiding model reduces to 
that of the ACDM model for n = oo. We find an analytic formula of the growth rate and 
the growth index using the quasi-static and sub-horizon approximation for the model with 
n 10. In section 5, we give a current constraint using the multipole power spectrum from 
the Sloan Digital Sky Survey (SDSS) luminous red galaxy (LGR) samples. A prospect of 
constraining the kinetic braiding model from a future redshift survey is also considered. In 
section 6, we investigate the behavior of the galileon field on small scales and the non-linear 
evolution of the density perturbation with the spherical collapse approximation. We also 
show that the kinetic braiding model with small n can be distinguished from the ACDM 
model using the galaxy cluster number count in a future survey. Throughout the paper, we 
use units in which the speed of light the Planck constant are unity, c = H = 1, and Mp\ is 
the reduced Planck mass related with Newton's gravitational constant by Mp\ = l/\/8irG. 
We adopt the Hubble constant Hq = 100/ikm/s/Mpc with h = 0.7, and the matter density 
parameter Vt^h 2 = 0.1344, unless stated otherwise. We follow the metric signature convention 
(-,+,+,+). 

2 Kinetic Braiding Model 

We consider the action that the galileon field <j) is minimally coupled to gravity, 



where R is the Ricci scalar, K((p, X) and G(4>, X) are arbitrary functions of the galileon field 
and the kinetic term X, which is defined by X = — g^S/ ^4>V v (J)/2, □</> = g^ v V ^Vy^ and £ m 
is the matter Lagrangian. The term G(4>,X)n\(f> is a general extension of the self-interaction 
term XO(f>, which is inspired from a decoupling limit in DGP model and the low-energy string 
effective action [64, 65]. We consider the extended model motivated by its phenomenological 
consequences, as is shown below. 

First, the evolution equation of the galileon field is the second order differential equation. 
The variation of the action with respect to g^ yields the equation for the gravity, 




in 



(2.1) 



(2.2) 



where T^ u is the energy momentum tensor of the matter and we defined 



- (V„GV^ + V M GV^) - G x n<j>V^V u <f>. 



(2.3) 



The galileon field equation is, 



Kj, + V a {K x V a 0) + 2G H X - 2G D0 
+ Gx4> (2XD(f> + 2 WV M V^V» 
+ G x [-(n^) 2 + (V M V^) 2 + iVV>V>] 
- G XX {VuXV u X - WV^V^VW) = 0, 



(2.4) 
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where we defined = dG((f>,X)/d4>, Gx = dG(c/>, X)/dX, and so on. In the present paper, 
we consider the case K((j), X) and G(<p, X) are the function only of X, as given by eqs. (2.14) 
and (2.15). In particular, in this case, the presence of the shift symmetry of the galileon field, 
cj) — > <p + c, allows the existence of a Noether current, which is given by 

J M = (K x - 2G - GxU&V^ + G X V M V^V>. (2.5) 

In this case, the galileon field equation (2.4) simply reduces to V M = 0. 

2.1 Background Evolution 

In this subsection, we summarize the background evolution of this model. For the spatially 
flat Friedmann-Robertson- Walker metric, 

ds 2 = -dt 2 + a 2 (t)S ij dx i dx j , (2.6) 

the gravity equations are given by 

3M^H 2 = p<j) + p m + Pr , (2.7) 
-M^(2H + 3# 2 ) =P4,+p r , (2.8) 

where H(= a/a) is the Hubble parameter, p m and p r are the energy density of matter and 
radiation, respectively, p r is the pressure of the radiation, and the energy density and the 
pressure of the galileon field are defined by 

P<f> = -K + K x tf - G^ 2 + 3G X H<£ 3 , (2.9) 
P<p = K - G^cj) 2 - Gx<i> 2 4>, (2.10) 

respectively. The galileon field equation is given by 

Kj, - (K x - 2G<p)(4> + 3H<P) - K^xtf - K X xW + G H 4> 2 

+ G X( p<p 2 (4> - %H4>) - 3G x {2H<fi4> + 3H 2 <p 2 + H<p 2 ) - 3G X xH<j) 3 4> = 0. (2.11) 

Note that the energy-momentum conservation for the matter and radiation remains the same 
as the standard one, 

Pm + 3H Pm = 0, (2.12) 
p r + 4:Hp r = 0. (2.13) 

In this paper, we consider the following case, 

K(X) = -X, (2.14) 

(2 \ n 

M 2 *) ' (2 ' 15) 

where n and r c are the model parameters, and r c is called the cross-over scale in the DGP 
model, which has a unit of length [71]. This model gives us phenomenologically interesting 
features, as we will show below. This model corresponds to the original galileon model in 
[46] if n = 1. 
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The Lagrangian is invariant under a constant shift symmetry (j> — > 4> + c, and the charge 
density of the Noether current (2.5) can be written as 



Jo = ^GxH - lj . (2.16) 

The galileon field satisfies Jo + 3HJq = and its solution is given by Jo oc 1/a 3 . Therefore, 
Jo approaches zero as the universe expands. It has been shown that the simplest attractor 
solutions is located at Jo = and has two branches, which are cj> = and 

(2.17) 



r 3G X H 

The case that = has ghostly perturbation as we will see in the next section and a self- 
accelerating solution does not exist. Therefore, we pick the other branch (2.17) throughout 
this paper. Then, we have 

G xx =2(n-l)-±j-, (2.18) 

Gxxx = 4(n - l)(n - 2)-i-, (2.19) 

SHcp 

and 

■■ 1 4>H 



2n-l H 



(2.20) 



Using (2.17) and the Hubble parameter as the present epoch (Hubble constant) Hq = H(ao) 
the cross-over scale is given by 



2n-l\ V 2 " 



3/? 



1 



6(1 - n - n 



rOj 



(2n-l)/4n 

Ho\ (2.21) 



where and £l r o are the density parameter of the matter and radiation, respectively, at 
present. In this case, we obtain the modified Friedmann equation, 



H \ 2 , „ x / H 



,i-n -n r0 ) — + Q a +O-0« ■ (2.22) 

-"0/ \-"0/ 
We also find that the effective equation of state of the galileon field can be written as 

WcS = J, = 3(2n - ijffl " 1 ( 2 - 23 ) 
6n + ^(a) (224) 



3(2n — fi m (a) — £l r (a)) 



where fi^(a) = p^/3M|j/J 2 , fJ m (a) = p m /3Afpj/J 2 , and O r (a) = pr/SM^H 2 , which satisfy 
O^(a) + fi m (a) + O r (a) = 1. In the derivation of the expression for w e s, we used 

H_ (2n-l)(3n m (a)+4n r (q)) 

/J 2 2(2n - fi m (a) - fi r (a)) ' 1 ' j 
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The effective equation of state has the asymptotic values, 



-— - — — — r (radiation dominated era), 
3 ( 2 ?7< 1 j 

^cfr — { (matter dominated era), (2.26) 

3 f 2?7/ 1 j 

— 1 (scalar field dominated era). 

The modified Priedmann equation (2.22) is known as the Dvali- Turner model, which 
was proposed as a phenomenologically modified model of the DGP model [67, 68]. At early 
stage of the cosmological expansion after the matter-dominated era, the matter energy density 
dominates the right-hand-side of the modified Friedmann equation (2.7), which reduces to the 
usual one, SM^H 2 ~ p rn . At late time, the matter energy density can be neglected, therefore 
the first term dominates the right-hand-side of the Friedmann equation (2.7). This term acts 
like the cosmological constant and asymptotically approaches the de Sitter expansion. This 
can be seen from that the effective equation of state approaches —1 for any n in the future. 

There are two special cases n = 1 and n = oo in background evolution. For n = 1, we 
can solve the Friedmann equation for H 2 , 



H\ 2 1 



H n J 2 



fl Q a~ 3 + f} r0 a" 4 + V (O a~ 3 + fi r0 a~ 4 ) 2 + 4(1 - fl - n r0 ) 



(2.27) 



Note that this solution is the same as the tracker solution in [43]. Here the matter density 
parameter is given by 

n m (a) = . 2 ^° Q " 3 ■ (2.28) 



n a~ 3 + ft r0 a~ 4 + y (^ a~ 3 + ^roa" 4 ) + 4(1 - Q - O r0 ) 

During matter-dominated era, fi m (a) — 1 gives w e g = —2. For large n, the Friedmann 
equation can be written as 



Hn 



~ 1 - fi - ^rO + ^oa" 3 + ^roo" 4 - (2.29) 



This is almost the same as the Friedmann equation of the ACDM model. 

The numerical calculation of the effective equation of state w c s for different parameters 
n is shown in figure 1 (left panel). As is expected, w e fr approaches —1 for large n. The right 
panel of figure 1 plots the comoving distance x( z ) = Jq dz' / 'H(z'). As n becomes large, the 
background evolution approaches that of the ACDM model. 

2.2 Stability Condition 

Let us study the stability of the model by perturbatively expanding the action of the galileon 
field at second order as 4> 4> + 6<f>. We neglect the metric perturbation, then the quadratic 
action becomes 



5S^ = I / d 4 Xy 



K (a)^ 2 -M( W) 2 



a- 



(2.30) 
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Figure 1. Left panel: The effective equation of state w e g as a function of redshift for ACDM (solid 
curve) and the kinetic braiding mode with n = 1 (dashed curve), n = 2 (dash-dotted curve), and 
n = 5 (dotted curve), respectively. Right panel: The comoving distance x(z), normalized by Hq, as 
a function of redshift for ACDM and this model. 



where 

«(o) = K x + K xx tf - 2G^ - G x ^ 2 

+ QG X H4> + 3G X xH<j) 3 + 12vrGG^ 4 , (2.31) 
/3(a) =K X - 2G <t> + G x ^ 2 + 2G X (4> + 2H4>) 

+ G X x<f> 2 4> ~ ^GG\^. (2.32) 

To avoid ghost and instability, we require k(o) > and c 2 s = j3{a)/n{a) > 0, respectively. 
One of the attractor solution, <fr = 0, has obviously ghostly perturbation since k = — 1 < 
in our model. Using the attractor condition (2.17), k(o), /3(a), and c 2 (a) can be written in 
terms of fi m (a) and fi r (a) as 

k(o) = (2n - fi m (a) - n r (a)), (2.33) 

_ n(50 m (a) + 6fi r (a)) - (fl>) + O r (a)) 2 

Pl j " 3(2n - n m (a) - O r (a)) ' 1 j 

2 _ ra(50 m (a) + 6» r (a)) - (n m (a) + O r (a)) 2 , 

CsW " 3(2n-fi m (a)-fi r (a)) 2 ' 1 ' 

Throughout the evolution of the universe, the matter density parameter O m (a) is always less 
than unity, therefore, we require n > 1/2 to avoid ghost instability, k > 0. In figure 2, the 
evolution of the sound speed of perturbations is plotted. The asymptotic behavior of the 
sound speed is given by 

{— - ; -7T (radiation — dominated era) , 
^In- 1 / (2-36) 
— -k (matter — dominated era) . 
3(2n — iy 

It is worthy to note the fact that the sound speed of the galileon field perturbation 
becomes zero for n = oo. In section 4, we show that the linear cosmological perturbation 
reduces to that of the ACDM model. In a case of large but finite value of n, the time- 
dependence of the galileon field perturbation will become important even on sub-horizon 
scales. For a case of small value of n, e.g., n 10, we can use the quasi-static approximation 
safely on small scales, which is demonstrated in section 4. 
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Figure 2. The sound speed of perturbations for n = 1 (dashed curve), n — 2 (dash-dotted curve), 
and n — 5 (dotted curve), from the top to bottom, as a function of scale factor. 



3 Constraint from Supernovae and CMB distance Observation 

The SCP Union2 Compilation is a collection of 557 type la supernovae data whose range of 
the redshift is 0.015 < z < 1.4 [72]. The distance modulus, which is the difference between 
apparent and absolute magnitude of the object, is given by 

"= 5iog (H) +25 - <"> 

where cIl(z) is the luminosity distance, 

d L {z) = (1 + z) f H-\z')dz' = (1 + z)x(z). (3.2) 
Jo 

The best-fit value of the present matter density Qq and n can be determined by x|n defined 
by 



XSN — 

i=l 



(p(Zj) - H6bs{ Zj)f 



(3.3) 



where Zi, fi bs( z i)i an d & bs( z i) are the redshift, the distance modulus, and the error of the 
i-th observed la SN. The left panel of figure 3 shows the confidence contour of A%sn on 
the plane fio and n for the kinetic braiding model. We find the minimum of the Xsn * s 
respectively given by Xsn = 542.7 for ACDM model and Xsn = 543.3 for the kinetic braiding 
model with n = 1. The kinetic braiding model with n = 1 requires the higher matter density 
parameter. A similar result is found in a recent paper [43]. 

The Wilkinson microwave background anisotropy (WMAP) observation also provide a 
constraint that relies on the distance to the last scattering surface. We here adopt the method 
using the shift parameter R, the acoustic scale I a, and the redshift of the decoupling epoch 
z* [73]. In the spatially flat universe, the acoustic scale and the shift parameter are written 
I A = nx{ z *)/ r s{z*) and R = sj 'Q,qHqx{ z *) > respectively, where x( z ) is the comoving distance 
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Figure 3. The left panel is the contour of Ax§ N on the plane fio and n for the kinetic braiding 
model. The dashed curve and the solid curve are the 1 a and 2 a confidence contours, respectively. 
The right panel is the same but of Ax^ MB . 



and r s is the comoving sound horizon at the decoupling epoch [74]. Then, we define the chi 
squared (see [4, 43] for details) 

where X{ = (Ja, R,z*), d\ = (302.09, 1.725, 1091.3), and the inverse covariance matrix 

/ 2.305 29.698 - 1.333 \ 

Cr: 1 = 29.698 6825.27 - 113.18 . (3.5) 
\ - 1.333 - 113.18 3.414 / 

In this analysis, we adopted Q r Q = 4.17 x 10~ 5 /i~ 2 . The right panel of figure 3 shows 1 a 
and 2 a confidence contours of Ax^ MB . Our result is consistent with [43] for the case n = 1. 
The higher n models better match the CMB distance observation. 

4 Linear Cosmological Perturbations 

In this section, we consider the linear evolution of cosmological perturbations. Hereafter, 
we consider the era after the matter domination. The metric perturbation in the conformal 
Newtonian gauge is given by 

ds 2 = -(1 + 2^)dt 2 + a 2 (i)(l + 2$)5 ij dx i dx j . (4.1) 

We write the perturbation of the galileon field by (f>(x, t) = <j)(t) + Scj)(x, t). Hereafter, <fi means 
4>{t). In appendix A, we summarized the perturbation equations in the general case. In our 
model, (2.14) and (2.15), the following equations are obtained. 
The (0, 0) component of the Einstein equation yields 

1 



2Ml x 



-3iJ($ - HV) + ^V 2 $ 



-K X 5X - G x ^30 3 6 - 12H<j) 3 ^ + 9H<}) 2 54> - ^ V 2 5(/>j 
3G X xH<j) 3 5X - dp, (4.2) 
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where we defined 



SX 



(4.3) 



and 5p is the matter density perturbation. The component of the Einstein equation 
yields 



2M 



pi 



(3H 2 + 2#)* + - $ - 3F$ 



+ 4</> 2 <^ 



and 



(4.4) 



(4.5) 



* + $ = 0. 
The (0, i) component yields 

2M 2 X - = -K X <PH - G x <j) 2 - H + 3H5(/)) + 5q, (4.6) 



where Sq describes the velocity field (see appendix A). The perturbed equation for the galileon 
field becomes 



-K x 



1 



30$ 2((j) + 3Hcf))y + 5(j) + 3HS(f> - -^V 



-SGxxxH^SX - G x 



12{ (H + 3H 2 )<}> 2 + 2H(t)(j) \V - + QH(j>5(t> 



3(j) 2 $ + 6(0 + 3H(t))(t)<S> - 9Hcl) 2 y 



-6^H4>+{H + 3H 2 
-Gxx 



+ 2H<f>)V 2 



3H(t)^ - 3{ 8H(f> 6 </> + (H + 3if 2 )0 4 }V 

12, 



SH^dcj) + 3<^ hH<F4> +(H + 3H 



-V 2 6(f> 



0. 



(4.7) 



The perturbed equation for the matter component is described in appendix A. Instead of 
solving the evolution equation of the matter perturbation, we use the following equation [34] , 
which is obtained by combining eqs. (4.2) and (4.6), 



+G X [-3$$ + 9H(j> 3 y - GH<p 2 5(j> - 9H 2 <p 2 5<j) 
-3GxxH<p 3 5X, 
where A c is the gauge invariant matter density contrast. 



(4. 
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4.1 Limit of n = oo 

In section 3, we showed that the background expansion of the kinetic braiding model ap- 
proaches that of the ACDM model as n becomes large. In this subsection, we show that the 
linear cosmological perturbation of the kinetic braiding model also reduces to that of the 
cosmological constant model for re = oo. 

From the attractor condition (2.17), in the limit of n = oo, we find that <f> is constant, 

0= [6(l-n )] 1/2 H M Ph (4.9) 

In this case, the leading term expanded eq. (4.7) in terms of 1/re gives 

5X + 3H5X = 0, (4.10) 

where we used (4.5) and 5X = 0<50— cfp"^ . Then, the solution of (4.10) is 5X = C/a 3 , where C 
is an integration constant, which is determined by initial condition of \l/ and 5(f). This solution 
is the decaying solution, therefore, one can set 5X = 0. In the case 5X = 0, the galileon field 
perturbations in the right hand sides of eqs. (4.4) and (4.6) disappear, and the equations 
reduce to those of the ACDM model. The evolution equations of the matter perturbations 
are not altered unless the gravitational potentials are altered. We are considering the initial 
conditions that is the same as the ACDM model. These suggest that the linear perturbations 
of the kinetic gravity braiding model for n = oo are equivalent to those of the ACDM model. 
In our numerical computation, we have confirmed that the cosmological perturbations of 
the kinetic braiding model with large re approach the result of the ACDM model (see also 
subsection 4.2). 

This property that our model reduces to the cosmological constant model can be un- 
derstood in the Lagrangian (2.1) with eqs. (2.14) and (2.15). In the case, (p = constant and 
5X = 0, K{X) and G(X) in the Lagrangian are constants, then the term GU(j> can be writ- 
ten as the total derivative, which can be dropped effectively when deriving the equation of 
motion. The term K(X) plays a role of the cosmological constant, and the kinetic braiding 
model reduces to the cosmological constant model for n = oo. This feature could be related 
to the fact that the sound speed of perturbations of the galileon field for re = oo is zero over 
the whole history of the universe. The perturbations of the galileon field can not propagate 
and does not induce an additional effect on gravity. 

The feature that the linear perturbations of the kinetic braiding model approaches those 
of the cosmological constant model appears for very large re, e.g., n ^ 1000, as will be shown 
in the next subsection. Such a behavior of large re model is different from that of smaller re 
model, e.g., n £ 10. For the case re £ 10, the quasi-static approximation can be adopted for 
the galileon field perturbations on the subhorizon scales, while it is not valid for re ^ 1000. 
This comes from the behavior of the sound speed of the galileon field perturbations. In the 
k-essence model with a zero sound speed, the scalar field perturbations cluster, and enhance 
the cosmological perturbations [75] , which is a contrast to our kinetic braiding model on the 
attractor solution (2.17). 

4.2 Numerical Calculation 

In this subsection, we show the evolution of the perturbations solved numerically. In figure 4, 
we plotted the evolutions of the gravitational potential (upper left panel), the perturbations 
of the galileon field 5(f) (upper right panel), and the density perturbations A c (lower panel), 
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for the kinetic braiding model with n = 1 with adopting different wavenumbers labeled there 
and for the ACDM model. Figs. 5, 6 and 7 are the same as figure 4 but with n = 10, n = 100 
and 1000, respectively. We can see that the perturbations approach those of the ACDM 
model as n becomes large, as we discussed in subsection 4.1. For n ^ 1000, the perturbations 
^ and A c with k ^, 0.01/iMpc" 1 are almost the same as those of the ACDM model, as one 
can see from figure 7. Since the sound speed of the perturbed galileon field become smaller 
as n increases, longer wavelength perturbations with larger n are not affected by the galileon 
field. Thus, the longer wavelength perturbations approach the ACDM model quickly as n 
increases. 

Let us discuss about the choice of the initial condition of the numerical computation. 
For the galileon field perturbations, the above results adopted the initial conditions 5(pi = 
and <50j = 0, where the subscript i denotes the initial values at ai = 10 -3 . In the Galileon 
scalar-tensor models [34], it has been pointed out that a different choice of the initial condition 
might give rise to an oscillating behavior in ^f, S(j), and A c . Figs. 8, 9 and 10 compares the 
gravitational potential vf, the perturbation of the galileon field 5(f>, and the growth factor 
A c , as a function of the scale factor, for the kinetic braiding model n = 1 with the different 
initial conditions 5(pi/M P \ = 0, Sfa/Mpi = 0.1^, Sfa/Mpi = ^ and 5<pi/Mp\ = 10^, 
respectively, where the other parameters are the same as those of Fig. 4. The behavior of 
the cosmological perturbations ^ and A c are almost same for k ^ 0.001/iMpc" 1 as long as 
\5(pi\/Mpi ^, 'Fj. The difference becomes smaller as the initial value \S(pi\/Mpi is smaller. 
Even for the cases 5<fii/Mp\ = and 5(f)i/Mp\ = lO^j, the the difference appears only for 
the large scales k ^ 0.0003/iMpc" 1 . We also checked that the different choice of 5^ does not 
alter this conclusion as long as {Scp^/HiMpi £ Vl/j, where Hi is the Hubble parameter at the 
initial time Oj = 10~ 3 . Therefore, the result is not very sensitive to the choice of the initial 
condition within a suitable range of it. 

If we assume that the galileon field undergoes an inflationary period at the early stage 
of the universe, fluctuations of order dfa ~ H- m f/2ir could be set, where H m { is the Hubble 
scale determined by the energy scale of the inflation. The condition \5(f)i\/Mp\ ^ yields 
Hmf/ '(2-7rMpi) ^, ~ 10 -5 , which might be satisfied as long as the energy scale of the 
inflation is not so high. 

For the model with larger n, the quasi-static approximation for the galileon field per- 
turbation fails. We can estimate the condition for the validity of the quasi-static approxima- 
tion, as follows. From the equation of motion obtained from (2.30), the condition is given by 
kS4> <C [3k 2 5(p/a 2 . Using the relations 5cp ~ H5(j) and 5cj) ~ H 2 5cp, the quasi-static approx- 
imation is valid as long as k 2 c 2 /a 2 3> H 2 . We find that the quasi-static approximation is 
valid for the wavenumber k £ 0.01/iMpc _1 for n ^ 10 from numerical results. However, the 
quasi-static approximation no longer works for n ^ 100. This is because the sound speed of 
perturbations of the galileon field approaches zero for very large n. 



4.3 Quasi-Static Approximation 

Hereafter, we consider the model with n ^ 10, in which small-scale perturbation equations 
can be expressed in simple forms by using the quasi-static and the sub-horizon approximation. 
Neglecting the time-derivative terms of perturbations and assuming that the wavelength of 
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Figure 4. The gravitational potential ^> (upper left panel), the perturbation of the galileon filed 
(upper right panel), and the growth factor divided by scale factor (lower panel) as a function of 
rcdshift for the ACDM model (thick dotted curve) for the kinetic braiding model of n = 1 with k = 
0.0003/iMpc" 1 , 0.001/iMpc" 1 , 0.005/iMpc" 1 , 0.01/iMpc" 1 , 0.05/iMpc" 1 , and the case of the quasi- 
static approximation (see subsection 4.3), respectively. The gravitational potential is normalized 
to unity at the initial time. 
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Figure 5. Same figure as figure 4 but with n = 10. 
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Figure 6. Same figure as figure 4 but with n = 100. 
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Figure 7. Same figure as figure 4 but with n = 1000. 



perturbation is well inside the Hubble horizon, 0(k cj/a ) 3> 0(H ), the evolution equation 
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Figure 8. The gravitational potential "J (upper left panel), the perturbation of the galileon filed 
(upper right panel), and the growth factor divided by scale factor (lower panel) as a function of 
redshift for the kinetic braiding model of n = 1 for k = 0.0003/iMpc" 1 , 0.001/iMpc^ 1 , 0.01/iMpc" 1 
with different initial conditions 6(f>i/Mp\ = 0.1 ^ (solid curve) and S(f>i/Mp\ = (dashed curve), 
respectively. We adopt the initial condition 8cf) i = 0. 



for the matter overdensity 5 in linear theory reduces to 

. V 2 
5 + 2H5- -=■*, 



(4.11) 



where we used the fact 5 = 5p/p ~ A c (see also appendix A). From eqs. (4.4) and (4.5), we 
have 



V 



T ^ ~ -47rGp<5 + AttGGx^ (f>—Tf, 

a z a z 



where we introduced tp(x,t) = 5(j>(x,t)/cj)(t). From eq. (4.7), we have 

p(a)—~-(p ~ -47rG — - — pd, 



(4.12) 



(4.13) 



for the galileon field perturbation. Combining (4.12) and (4.13), we obtain the modified 
Poisson equation, 



VP ~ 4irGesp5, 



(4.14) 
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where the effective gravitational constant is given by 

G e s = G 



1 + 4, G G ^ 4 



/3(a) 

_ 2n + 3nQ m (a) - £l m (a) . . 

n m (a)(5n-n m (a)) ' 1 ' j 

Here, in the second line, we used the attractor condition (2.17). Then, the evolution equation 
for the matter overdensity in linear theory can be written as 

5 + 2H5 ~ 4,nG cSP 5. (4.16) 

The crucial difference between the kinetic braiding model and the scalar-tensor galileon the- 
ory [33-35] is the absence of an effective anisotropic stress in the right-hand-side of eq. (4.5). 
The effective gravitational constant G e g is nearly equal to G at early times and becomes larger 
than G at late times. (See the dashed curve in the right panel of figure 17). The enhancement 
of the effective gravitational constant leads to an enhancement of the growth factor of matter 
density perturbations. Although the background evolution for large n approaches the ACDM 
model, the growth history is different due to the time-dependent effective gravitational con- 
stant. The left panel of figure 11 compares the evolution of the growth factor divided by the 
scale factor, 5 /a, which is normalized as 1 at early stage of the evolution. The growth factor 
of the kinetic braiding model is larger than that of the ACDM model throughout its growth 
history. In the right panel of figure 11, the linear growth rate f(a) = din 5 / '(fin a is plotted. 
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Same figure as figure 
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Figure 11. Left panel: The growth factor divided by scale factor as a function of scale factor for the 
ACDM model (solid curve) and the kinetic braiding model n = 1 (dashed curve), n = 2 (dash-dotted 
curve), and n = 5 (dotted curve), respectively. Right panel: The linear growth rate as a function of 
scale factor. 



4.4 Analytic Formula for Growth Index 

The growth rate can be parametrized by the growth index 7(a) defined by f(a) = f2 m (a) 7 (°). 
It has been known that the growth rate of the ACDM model is well approximated by the 
nearly constant value 7 = 6/11. This parametrization was first introduced in [76] and has 
the potential to distinguish dark energy models and modified gravity models [77]. In this 
section, we derive the analytic formula of the growth index for the kinetic braiding model. 
As we will show below, the behavior of the growth index is unique in the kinetic braiding 
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model, which will be important to distinguish between the kinetic braiding model and the 
ACDM model when n is not infinite. 

The matter evolution equation (4.16) can be written in terms of the growth rate /(a) = 
din 5 /din a as 

+ / 2 +(2 + 4V^ m (a) = 0. (4.17) 



dlna \ H 2 J 2 G 

By using the background equation and (2.17) , we have 

H 3(2n - l)Sl m (a) 

IP ~ ~ 2{2n - fi m (o)) ' 
d£l m (a) 6nSl m (a)(l — Sl m (a)) 
din a 2n — Sl m (a) 

Using eqs. (4.15), (4.18) and (4.19), eq. (4.17) can be written as 

6nfi m (a)(l - fi m (a)) f dj(a) 7(a) 



(4.18) 
(4.19) 



2n — Sl m (a) 



{^Llnn m (a) + ^L} +fimW ^) 
[dSZ m (a) ii m (a) J 



8ra-6nSl m (a) -S7 m (a) (a) _ 3(2n + 3nSl m (a) - Q m (a)) 
+ 2(2n-O m (a)) mi j 2(5 - fi m (a)) " 1 j 

We assume that the growth index can be approximated in the following expanded form, 

m 

7(a) =<E 0(1 (4-21) 

Then, we expand eq. (4.20) around Sl m (a) = 1 and we have the asymptotic value of the 
growth index at high redshift, 

Co = 3n l 16n " 5) ■ (4.22) 
s llOra 2 - 47n + 5 y ' 

The higher order coefficients can be easily found. By way of example, table 1 summarizes the 
coefficients for the model n = 1, 2, 3, 4, and 5. The redshift evolution of 7(a) is important 
in the kinetic braiding model. This is shown in figure 12, which plots the growth index 7(a), 
by solving eq. (4.17) numerically, as a function of redshift, for n = 1, 2, and 5, respectively. 
Figure 13 compares the numerical result and the approximate formula of (4.21) of the model 
n = 1 as a function of redshift. The left panel of figure 13 plots the numerical result and the 
approximate formula with m = and m = 5, respectively, where m is the maximum index 
in (4.21). The right panel shows the relative error (/ — Sl m (a) 7 ^ a ^)// for m = 0, 1, 2, 3, 4, 
and 5. The approximate formula (4.21) with m = 5 reproduces the exact growth rate / to 
better than 0.46%, and the growth factor can be obtained in a highly accurate manner, by 
D\(a) = aexp [J Q a din a' {Sl m (a') 7 — 1}]. For the other case in table 1, we have the similar 
result and confirmed that the growth factor can be well approximated by the formula up to 
m = 5. 

On very small scales, the non-linear terms of (p become more important than metric 
perturbations in the case of small n. The field equation for the galileon field including the 
non-linear terms is given by 

P{a)%y + ^ [(V,V^)(V l VV) - (VV) 2 ] = -inG^fpd, (4.23) 
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Figure 12. The growth index of the kinetic braiding model as a function of redshift for n = 1 (dashed 
curve), n = 2 (dash-dotted curve), and n = 5 (dotted curve), from the top to bottom, respectively. 
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Figure 13. Comparison between the numerical result and the approximate formula of (4.21) for the 
model n = 1. Left panel is the growth index as a function of redshift. The solid curve is the exact 
result from numerically solving the evolution equation, but the other two curves are from the formula 
in the expansion (4.21) with m = 1 (dotted curve) and m = 5 (dashed curve), respectively. Right 
panel is the relative error Sf/f for the case when the maximum value m is set as 0, 1, 2, 3, 4, and 
5, respectively. 





n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


Co 


0.48529 


0.46154 


0.45316 


0.44895 


0.44643 


Ci 


-0.03373 


-0.03153 


-0.02922 


-0.02791 


-0.02708 


c 2 


-0.02814 


-0.02037 


-0.01876 


-0.01805 


-0.01764 


C3 


-0.01667 


-0.01470 


-0.01393 


-0.01355 


-0.01333 


a 


-0.01367 


-0.01174 


-0.01119 


-0.01093 


-0.01078 




-0.01106 


-0.00976 


-0.00936 


-0.00916 


-0.00905 



Table 1. Coefficients Q for the model of n = 1, 2, 3, 4, and 5. 



which holds for the general Lagrangian (2.1) as long as the quasi-static approximation works. 
The equations for the metric perturbation remain the same as those of the linear perturbation 
theory. From eq. (4.23), the non- linear terms in the galileon field equation become dominant 
when ip > fia 2 /k 2 Gx4>- From eq. (4.13), we may estimate (p ~ a 2 H 2 Gx<p 2 S/k 2 (3(j), by using 
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the relation, H 2 ~ p m /3Mp, . Combining these relations, the non-linear terms become 
dominant when 5 > ft 2 / H 2 G 2 X 4> 2 ~ 0(1). Therefore, for 6 > 1, we need to take into account 
the non- linear terms, and these terms play an important role in the Vainshtein mechanism 
as we will see in section 6. 

In the limit of n = oo, the linear cosmological perturbation approaches the ACDM 
model. It would be interesting to clarify what is the nature of the gravity in the model 
with large n on the small scales where the Vainshtein mechanism works in the static limit. 
However, this is outside of the scope of the present paper, and we only consider the case 
n 10 in the latter part of this paper. 

5 Cosmological Constraints from the Redshift-space Distortion 
5.1 Current Constraint 

The multipole power spectrum is useful for measuring the redshift-space distortion, which 
plays a vital role of testing gravity on the scales of cosmology [78-83]. In this section, 
utilizing the monopole and quadrupole spectra from the SDSS LRG sample of the data 
release 7 [31, 84], we investigate a constraint on the kinetic braiding model. The multipole 
power spectrum is defined by 

P(k,n)= Pl(k)W(2l + l), (5.1) 

Z=0,2,4,... 

where P(k,/i) the redshift-space power spectrum, Ci(n) are the Legendre polynomials, fi is 
the directional cosine between the line of sight direction and the wavenumber vector. The 
monopole Po(k) represents the angular averaged power spectrum, and Piik) is the quadrupole 
spectrum which gives the leading anisotropic contribution. For the theoretical modeling of 
the redshift-space power spectrum, we adopt the fitting formulas developed by Jennings et 
al. [85] and Peacock and Dodds [86], with the transfer function by Eisenstein and Hu [87]. 
We consider the galaxy redshift-space power spectrum [88], 

P ga x(k,fi) = (b 2 (k)P ss (k) + 2fb(k)P &e {k)^ 2 + f 2 P ee (k)^) e~^ k ^ 2 , (5.2) 

where Pss(k) is the nonlinear matter power spectrum, Pgg^k) is the power spectrum of the 
velocity divergence, and Psg(k) is the cross power spectrum of matter and the velocity diver- 
gence, b(k) is the clustering bias, and a v is the velocity dispersion. Jennings et al. proposed 
a fitting formula for the redshift-space power spectrum of the form (5.2), assuming b{k) = 1. 
The fitting formula relates the nonlinear matter power spectrum P$${k) to P$e(k) and Pgg(k). 
Although the accuracy of the fitting formula for the kinetic braiding modified model is not 
verified, we assume its validity, and use it in this section. Also we consider the scale dependent 
bias in the form, 

6(t) = 6 « + 6l (al^) fe ' (53) 

where bo, b±, and 62 are the free parameters. The multipole power spectrum is measured 
assuming the fiducial distance-redshift relation s = s(z) of the spatially flat ACDM model 
with f2o=0.28. In order to compare the theoretical prediction of the kinetic braiding model 
with the observed power spectrum, we must take the cosmological redshift-space distortion 
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into account [89] . Denoting the comoving distance and redshift relation of the kinetic braiding 
model, x = x( z )i the galaxy power spectrum with the cosmological redshift-space distortion 
is 

w = j_ P ^ k ^ + _ m y ( , 4) 

where we defined cm = dx(z)/ds(z), c± = x( z )/ ' s ( z )i an d we use the mean redshift z = 0.3 of 
the SDSS LRG sample. 

We use the monopole and quadrupole power spectra in the range of wavenumber, 
0.02/iMpc^ 1 < ki < 0.2/tMpc~ , and compute the chi squared, 

Xrd = E E m^-Pt^k^C^iKkjXPeik^-P^ikj)), (5.5) 

i,j £,£'=0,2 

where Pf° s (ki) is the observed power spectrum and C^{ki,kj) is the covariance matrix. 
Pf as {ki) and Cm (hi) are from refs. [84], which used luminous red galaxy sample of the Sloan 
Digital Sky survey data release 7 with the contiguous part of 7150 deg 2 sky coverage in the 
redshift range 0.16 < z < 0.47, and the corresponding 1000 mock catalogues. 

Figure 14 shows the contours Ax^ D on the Qo — °8 plane for the ACDM model (upper 
left panel) and for the kinetic braiding model with n = 1 (upper right panel) , n = 2 (lower left 
panel) and n = 5 (lower right panel), respectively. In each panel, the la and 2a confidence 
contour-levels are given by dashed curve and solid curve, respectively. In the analysis, we 
fixed the spectral index n s = 0.96 and the density parameter for baryon £l\yh 2 = 0.0225, 
and the Xrd ^ s computed to minimize (5.5) by fitting the bias parameters bo, b±, hi and 
the velocity dispersion a v . In this figure, we may understand that the constraint on fio 
comes from the baryon acoustic oscillation feature, and the constraint on as comes from 
the amplitude of the monopole and the quadrupole spectrum. From figure 14, the kinetic 
braiding model with smaller n favors the smaller values of f^o and as, compared with the 
ACDM model. However, the difference is not quite large, and the constraint is not very tight. 

Our result is obtained by extrapolating the quasi-nonlinear formula to the kinetic braid- 
ing model. The Vainshtein mechanism is important on small scales in the kinetic braiding 
model, as in the case of the DGP model [90-92]. In the present paper, however, we ne- 
glect the nonlinear Vainshtein effect because we need to consider only rather large scales, 
k 0.2/iMpc" 1 . But, the validity will be necessary to be checked using N-body simulations 
in future. 

5.2 Future Prospect 

We now consider a future prospect constraining the kinetic braiding model using the multipole 
power spectra. Let us consider a redshift survey like the SuMIRe (SUbaru Measurement of 
Imaging and REdshift of the universe) PFS (Prime Focus Spectrograph) survey [93] , which 
assumes the survey parameters like those of the WFMOS survey [94]. In the SuMIRe PFS 
survey, we assume the sky coverage is 2000 deg 2 over redshift range 0.8 < z < 1.6 with the 
averaged number density of galaxy n = 3 x 10~ 4 /i 3 Mpc~ 3 . 

Figure 15 shows the monopole (upper panels) and quadrupole (lower panels) power 
spectra multiplied by k, kPo(k) and kP2(k) for the fiducial ACDM model with Q^h 2 = 0.1344, 
and for the kinetic braiding model with Qq = 0.27, 0.28, and 0.29. Here, we fixed n = 5. The 



- 21 - 



0.9 



0.8 



b 0.7 
0.6 
0.5 

0.9 
0.8 
b 0.7 
0.6 
0.5 







/ y 
/ / 


\ 


/ / 


\ 


/ / 




i 1 
( 






; 


\ \ 


/ 


(a) ACDM 



0.26 0.28 0.30 0.32 0.34 
fi 



(c) n=2 




0.26 0.28 0.30 0.32 0.34 




0.9 
0.8 
b 0.7 
0.6 
0.5 



0.26 0.28 0.30 0.32 0.34 



(d) n=5 




0.26 0.28 0.30 0.32 0.34 
O 



Figure 14. Contour of ^-Xrd on the — °8 plane. The upper left panel (a) is for the ACDM model. 
The upper right panel (b) is the kinetic braiding model with n = 1, the lower left panel (c) is n — 2, 
and the lower right panel (d) is n = 5, respectively. In each panel, the dashed curve and solid curve 
are the ler and 2a confidence contour-levels, respectively. 



left panels assume the subsample of the the range of redshift 0.8 < z < 1.2, while the right 
panels assume the subsample of the range redshift 1.2 < z < 1.6. The theoretical modeling 
of the power spectrum is the same as that in the above subsection, which is evaluated at the 
mean redshift of each subsample. In computing the theoretical power spectrum, we fixed the 
initial amplitude of the fluctuation so that the matter power spectrum of the ACDM model 
gives ag = 0.8. Other cosmological parameters are fixed as n s = 0.96 and Q^h 2 = 0.0225. The 
cosmological redshift-space distortion is taken in the power spectrum of the kinetic braiding 
model into account, assuming that the distance-redshift relation of the ACDM model is 
adopted in data analysis. The error bars for the ACDM model is evaluated by [95] 



(2vr) 3 
AT4 



dfi 



Vn 2 [l + nP(k,n)] 



-2 ' 



(5.6) 



where 14 is the volume of a shell in the Fourier space, and V is the survey volume, V = 
A r 2max dz(dsldz)s 2 , where A is the survey area, and z max and z m ; n are the maximum and 

J *min 

minimum redshifts of the survey, respectively. 
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Figure 15. Monopole (upper panels) and quadrupole (lower panels) power spectra multiplied by k, 
kPo(k) and kP2(k) for the fiducial ACDM model and for the kinetic braiding model of n = 5 with 
flo = 0.27, 0.28, and 0.29. The left panels are the subsample of the range of redshift 0.8 < z < 1.2, 
while the right panels are the subsample of the range redshift 1.2 < z < 1.6. The error bars is 
estimated using eq. (5.6). 

As for the bias b(k), We adopted the scale-independent bias b = 1.5 and the velocity 
dispersion a v = 350km/s for the ACDM model, but the scale-dependent bias (5.3) for the 
kinetic braiding model. In computing the power spectrum of kinetic braiding model, we 
assumed bo, b\, b<i and a v as free parameters which are determined so as to minimize the 
value 

E E l P i ACBU ( k i) - Pf n ^° s (k i )] 2 /AP i (k i ) 2 , (5-7) 

i 1=0,2 

where P e ACDM (h) and Pf raiding (^) are the power spectrum of the ACDM model and the 
kinetic braiding model, and we used 0.02/iMpc" 1 < h < 0.2/iMpc" 1 . 

Figure 15 shows that the monopole spectrum can be almost same in the different models 
because of the bias, but the quadrupole spectrum is different. Thus, the combination of the 
monopole and quadrupole power spectra is useful to distinguish between the ACDM model 
and the kinetic braiding model, depending on the model parameters. In 5.3, we estimates a 
constraint on SIq and n at a future quantitative level with the Fisher matrix analysis. 
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5.3 Fisher Matrix Analysis 

In this subsection, let us discuss the future prospect of a constraint with galaxy power 
spectrum quantitatively using Fisher matrix analysis. We adopt the following Fisher matrix, 
for simplicity, 

j_ /•*— y+i a(iffi(M))a(i& R (M)) y 
y " A min ^ «>« ^ (PCR (j fc, M ) + l/n)2' 

(5.8) 

where V is a survey volume, n is a mean number density of galaxies, Q{ denotes a model 
parameter. In this Fisher matrix analysis, we consider the 6 parameters n, Qq, bo, b\, 62 
and a v . The amplitude of the initial cosmological perturbation is fixed so that as = 0.8 for 
the ACDM model. We also fixed the spectral index n s = 0.96 and Q b h 2 = 0.0225. The 
curves in figure 16 show the 1-sigma (dashed curve) and 2-sigma (solid curve) confidence 
contours in the n and f^o plane, by marginalizing the Fisher matrix of the 6 parameters, n, 
^0) ^0; b\, 62 and cr v over the parameters other than n and Qq. In figure 16, we assumed 
the redshift survey like WFMOS/SuMIRe PFS survey, adopted in the previous subsection. 
The sky coverage is 2000 deg 2 , the minimum and the maximum redshifts are z m - m = 0.8 and 
z max = 1-6, respectively, and the averaged number density of galaxy is n = 3x 10~ 4 /i 3 Mpc -3 . 
The target model is the kinetic braiding model with n = 5, Qoh 2 = 0.1344, &o = 2, 61 = 0.1, 
62 = 0.1, and o v = 350km/s. We assumed that the power spectrum is measured adopting the 
fiducial distance-redshift relation s = s(z) of the spatially flat ACDM model with f2o=0-28. 

Figure 16 implies that the constraint on Qq is tight, while the constraint on n is not 
very tight. However, the 1-sigma (2-sigma) error in determining n is An ~ 5 (An ~ 10) 
when the target parameter is n = 5. Then, it is possible to distinguish between the kinetic 
braiding model and the ACDM model if n ^ 10. However, the constraint on n becomes weaker 
when the target value of n becomes large, since the background expansion and cosmological 
perturbations in the kinetic braiding model for n ^ 1000 is almost the same as the ACDM 
model. Then, it will be difficult to distinguish between the ACDM model and the kinetic 
braiding model by a measurement of the redshift-space power spectrum when n is large. 

6 Non-linear Density Perturbation 

In this section, we investigate the nonlinear effect of the perturbed galileon field on small 
scales. As discussed at the end of section 4, the non-linear term in the perturbed galileon 
field equation (4.23) becomes important for 5 ^ 1. Therefore, when considering the non- 
linear cosmological perturbations, we need to solve eq. (4.23). However, it is not easy to 
solve eq. (4.23) in general. For simplicity, we assume that the density perturbation 5p has a 
spherically symmetric distribution and a top-hat profile. Then, the perturbation equation for 
the galileon field can be solved analytically, and it gives a qualitative picture of the Vainshtein 
mechanism in the kinetic braiding model. Note that eq. (4.23) is identical to the equation 
for the brane bending mode <p in the DGP model [91, 96, 97] and the full profile of the brane 
bending mode (f is investigated in [98]. The behavior of <p in a top-hat density profile remains 
the same in the kinetic braiding model. 
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Figure 16. The 1-sigma (dashed curve) and 2-sigma (solid curve) contours in the flo — n plane when 
we assumed the WFMOS/SuMIRe PFS survey, A = 2000 deg 2 , the range of the redshift 0.8 < z < 1.6 
and the averaged number density of galaxy n = 3 x 10 _4 /i 3 Mpc~ 3 . The target model is the kinetic 
braiding model with n = 5, VLqIi 2 = 0.1344, bo = 2, b\ = 0.1, 62 = 0.1, and a v — 350km/s. We also 
fixed Vt^h 2 = 0.0225, n s — 0.96 and the amplitude of the initial cosmological perturbation so that 
cr 8 = 0.8 for the ACDM model. 




We start with recasting the Poisson equation (4.14) and the field equation (4.23) by 
replacing the comoving nabla operator with the physical nabla operator, then we have 

V 2 tp + A 2 (a) [(VjV.^)(V*VV) - (VV) 2 ] = -4^GC(a)5p, (6.1) 
V 2 ^ = 4ttG P 5 - £(a)vV (6.2) 

where we defined 



A 2 (a) 



G x <j> 



C(a) 



G x <t> 2 



/3(a)' * w f3(a)<p 
As we will show below, the following combinations are useful, 

(1 - n m (a))(2n - O m (a)) 



e(a)C(a) = 
A 2 (a)C(a) 



S7 m (a)(5ra - O m (a)) 
2n — Sl m (a) 



(6.3) 

(6.4) 
(6.5) 



H(a)£l m (a)(5n — f2 m (a)) 

Under the assumption that the density perturbation 5p(r) is spherically symmetric, the field 
equation (6.1) can be written as 



]_d_ 

r 2 dr 



,dip 
dr 



2A 2 (a) d 
r 2 dr 



dip 
dr 



-4nG((a)6p. 



(6.6) 



Then, we define the enclosed mass perturbation, m(r) = 4tt r l2 5p(r')dr' . Integrating 
eq. (6.6) over r 2 dr, we have 

r 2 ^-2A 2 (a)r(^ 
dr \ dr 



-Q{a)Gm{r). 



(6.7) 
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Solving for dip/dr, we have 



dip 
dr 



ry 



4A 2 (a) 



9 — 

ry 



(6. 



where g(x) 



and the Vainshtein radius of the enclosed mass is 

1/3 



r v = [8\ 2 (a)((a)Gm{r)] 



1/3 



8G 2 x <f> 2 



Gm(r) 



(6.9) 



Practically, eq. (6.7) yields two branches of solutions, which correspond to the attractive 
and repulsive force. We require that the solution matches the linear solution, given by 
dip/dr oc r~ 2 at r 3> ry, hence we pick the solution whose value vanishes at an infinite 
distance. The Vainshtein radius for a point source can be estimated using eq. (6.9) and 
(2.17), 



4r„ 



ry 



(~9# 2 /3 2 
1.91 x 10 2 



1/3 



2n — Qmta) 



£l m (a)(5n - O m (o)) 



2/3 



H 



-2/3 



M 



M 



1/3 



pc, 



(6.10) 



where r g is the Schwartzschild radius of the point mass. Using the present values of H and 
£l m for n = 1, the Vainshtein radius is approximated by ry = 2.3 x lO 2 (M/M0) 1 / 3 pc. The 
asymptotic behavior of the galileon field well inside and outside the Vainshtein radius is 



dip 
dr 



-C(«) 



Gm 



3/2 



4A 2 (a)rV2 



r » ry, 



r <C ry, 



(6.11) 



respectively. Well outside the Vainshtein radius, r 3> ry, dip/dr is proportional to the 
Newtonian potential and this solution corresponds to the linear solution. On the other hand, 
well inside the Vainshtein radius, r ry, the galileon field <p approaches nearly constant 
[98] and does not contribute to the modified Poisson equation (6.2). 

6.1 Top-Hat Profile 

Hereafter, we consider the following top-hat density profile p(r) or, equivalently, the mass 
perturbation m(r) 



p(r) 



p + 5p 
P 



r <R 
r>R' 



m(r) 



f 8M(r/Rf r<R 
\SM r>R 



(6.12) 



We neglect any compensation underdensity swept out by the prior evolution of the top hat, 
and assume that a density enhancement is constant. For r < R, dip/dr is proportional to r. 
Then, the non-linear term in eq. (6.1) can be written as 



(W^)(v { v^) - (vV) 5 



-^(vV) 2 . 



(6.13) 
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Therefore, eq. (6.1) can be written as 



2- 2 ^!M(VV) 2 = -47rGC(a)8p. 



(6.14) 



Setting r = R and solving for V 2 (/?, we have 



VV=-8,G««)(^I 9 I^Up, 



1/3 

where i?y = [8A 2 (a)C(a)G5Af] . Then, the modified Poisson equation (4.14) becomes 



(6.15) 



where the non-linear effective gravitational constant is given by 



47rG^ ff L = AttG 



1 + 2 »)l|) '(£ 



(6.16) 



(6.17) 



The first term of the right-hand-side of eq. (6.17) is the contribution from the Newtonian 
gravity and the second term induces the additional attractive force due to the perturbed 
galileon field. In the case R — > 0, corresponding to 5p —> oo, the effective gravitational 
constant approaches the standard gravitational constant G and this is the Vainshtein 
mechanism in the spherical collapse in the cosmological perturbations (See the solid curve in 
the right panel of figure 17). 

6.2 Spherical Collapse 

In this section, we examine the spherical collapse which is a powerful tool for understanding 
the formation of a bound system in non-linear regime [99-104]. The evolution of a spherical 
collapse undergoes three phases: (i) turn around, (ii) collapse, and (iii) virialization. A 
spherical overdense region expands along the expansion of the universe up to the maximum 
radius as the density perturbation evolves proportionally to the scale factor at early times. 
Then, it starts to shrink toward the center of the spherical region and finally collapses. The 
spherical collapse approach enables us to evaluate the critical density contrast 5 C and the 
virial density A v ; r , which need in predicting the mass function of the halos. We assume that 
perturbed matter distribution remains a top-hat profile and each shell does not cross during 
collapse. 

We start with the fully non-linear evolution equation for the matter over-density, which 
can be obtained from the energy-momentum conservation, 



4 8 2 



+ 2H5 = (1 + 5)V 2 ^. 



31 + 5 

We assume that the total mass inside R is conserved during collapse, 



M 



4:7T 4:7T 

— R 3 p = — R 3 p{^ + $)= constant. 
3 3 



Differentiating eq. (6.19) with respect to time gives 



■ 4 5 2 
6 + 2H5--, 



31 + 5 



(6.18) 



(6.19) 



(6.20) 
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Combining (6.18) and (6.19), the evolution equation for the spherical collapse is given by 



R 3 

= H 2 + H ^-pS, (6.21) 



as 

R 1 



where we used (6.17) in the second line. For convenience, we rewrite the equation of motion 

3 [V 2 * b + V 2 * g ] , (6.22) 

where we have defined the gravitational potential contributed from the effective background 
term V 2 ^b = —3(H 2 + H) and from the effective gravitational constant V 2 ^ g = AirG^pS. 

6.3 Virial Theorem 

For a spherically symmetric top-hat model, the total potential energy from the effective 
background term and the gravity is given by 

<"=S>. = -a"jf£rE£. (0.23) 

s=b,g s=b,g 

where the potential energy of the gravity and the effective background term are 

3 G^MSM 

Wg = "5^ ' (6 - 24) 

W h = l(H 2 + H)MR 2 , (6.25) 
5 

respectively. The virial theorem can be obtained by integrating the Boltzmann equation 
over all position, which gives 2T v ; r + W v - U = 0. The virial radius i? v ; r is usually calculated 
by requiring the energy conservation at turnaround and virialization. However, there is the 
energy non-conservation problem in the dark energy models, DGP and other modified gravity 
models [98, 105]. The violation of energy conservation also arises in the kinetic braiding model 
due to the time- varying and G e s. Hence, we need carefully choose the condition for the 
virialization epoch. In the present paper, we follow the approach in [98]. We choose the 
virialization epoch so that the virial condition satisfied, 2T v i r (a v i r ) + W(a v j r ) = 0, where the 
kinetic energy during collapse can be computed by 

T = f d z x\pv 2 = -^-MR 2 . (6.26) 
J Z 11) 

We define the virial radius i? v ir as the radius at the virialization epoch. The virial density is 
given by 

A vir = = [i + 5(R vil )} f^^l\ 3 , (6 .27) 

Pcollapse \ ^"vir / 

where p v - a = p(l + <5)| a=a . is the density at the virialization epoch and /? C oiiapse is the 
background matter density at the time of collapse. 
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We are now ready to solve the spherical collapse eq. (6.21), numerically. The total mass 
conservation within R gives 



\ 3 / o \ -3 
a \ I R 



^°»=UJ UJ (6 - 28) 

where i?; is the initial radius of the perturbation, 5i is the initial density fluctuation, and a\ 
is the initial scale factor. We set the initial scale factor a; = 10~ 5 . Turnaround and collapse 
are, respectively, defined by R = and R = 0. The left panel of figure 17 shows the evolution 
of the scaled radius for the ACDM and the kinetic braiding model with n = 1, 2 and 5. We 
set the initial density perturbation 5{ for the ACDM so that collapse occurs at a = 1. For 
the same initial density perturbation <5;, the spherical overdensity region collapses earlier in 
the kinetic braiding model compared with the ACDM model due to the additional attractive 
force. The right panel of figure 17 shows the linear and non-linear effective gravitational 
constant as a function of redshift. The radius of the overdensity region enters the Vainshtein 
radius around a ~ 0.48 for the case n = 1, and approaches the standard gravitational 
constant G in the kinetic braiding model. Thus, the non-linear effect restores the general 
relativity in high density regions through the Vainshtein mechanism. Figure 18 shows the 
evolution of the critical density contrast 5 C and the virial density as a function of redshift. The 
critical density contrast is the value of the linear-theory density contrast reached to infinite 
density when the nonlinear top-hat perturbation collapses. One can see that 5 C for the kinetic 
braiding model approaches that of the Einstein de Sitter universe, 5 C = 3(12vr) 2 / 3 /20, at high 
redshift. The virial density for the kinetic braiding model also approaches the value for 
Einstein de Sitter universe, A v ; r = 18-7T 2 , at high redshift. This is because the evolution of 
background and cosmological perturbations in the kinetic braiding model is identical to that 
in Einstein de Sitter universe at early times. Table 2 shows the numerical results of the 
critical density contrast 5 C and the virial density A v ; r at z = for the ACDM model and the 
kinetic braiding model with n = 1, 2, 3, 4, and 5. 





ACDM 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 




1.675 


1.696 


1.705 


1.708 


1.710 


1.711 


A v ir 


360.4 


305.8 


306.2 


306.3 


306.3 


306.4 



Table 2. 5 C and A v j r for the kinetic braiding model with n = 1, 2, 3, 4, and 5, respectively, at 
z = 0. 



6.4 Galaxy Cluster Counts 

We here give a simple prediction about the mass function of halos by using the Press-Schechter 
formalism [106]. The number density of virialized clusters in an unit mass range as a function 
of redshift is given by 



dn(M, z) 
dM 



2_p_ 

vrM 2 



S c (z) dlna(M,z) 
a(M, z) dlnM 





1 


exp 


2a 2 (M,z)_ 



(6.29) 



where M = 47ri? 3 p/3, and cr(M, z) is 



a 2 (M,z) = ^l k 2 dkP(k,z)W(kRf 



(6.30) 
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Figure 17. Left panel: Evolution of the scaled radius R/R[ of perturbations as a function of the 
scale factor for ACDM (solid curve) and the kinetic braiding model with n = 1 (dashed curve), n = 2 
(dash-dotted curve), and n = 5 (dotted curve), respectively. Right panel: Evolution of the effective 
gravitational constant G c s in the linear theory (dashed curve) of (4.15) and the non-linear theory 
(solid curve) of (6.17), normalized by G, as a function of the scale factor for the kinetic braiding 
model with n = 1 (dark black curve) and n = 5 (light red curves). 
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Figure 18. Left panel: The rcdshift evolution of the linear overdensity for ACDM model (solid 
curve), Einstein dc Sitter universe (thick-dashed line), and the kinetic braiding model as a function 
of redshift for n = 1 (dashed curve), n = 2 (dash-dotted curve), and n = 5 (dotted curve). Right 
panel: The redshift evolution of the virial density. 



with the top-hat smoothing function W(kR) = 3(sin kR — kRcos kR)/(kR) 3 and the linear 
matter power spectrum P{k, z) at the redshift z. The number count of halos in unit redshift 
within the range of mass, M\ < M < M2, is given by integrating the halo mass function, 

dN dV f Ah dn(M,z) 

J dM, (6.31) 



dz dz J Ml dM 

where dV/dz = A\ 2 {z)dx/dz is the comoving volume of the unit redshift and A is the survey 
area. 

We consider the expected redshift distributions of clusters in the upcoming survey, the 
Southern Polar Telescope (SPT) Sunyaev Zeldovich (SZ) survey with a limiting flux density 
of fuoiim = 5mJy at uq = 150GHz and a sky coverage of A = 4 x 10 3 deg 2 . The relation 
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Figure 19. The expected number of clusters N(M > Mr im ) in the redshift bins of Az = 0.1, in 
the SPT survey for the ACDM model with as = 0.8 and the kinetic braiding model n — 5. The left 
panel compares the ACDM model with as = 0.8 and the kinetic braiding model with as = 0.75, 0.8, 
and 0.85, respectively, with fixing HqH 2 = 0.1334 for all cases. The right panel compares the ACDM 
model with the kinetic braiding model with flo = 0.27, 0.28 and 0.29, respectively, with fixing the 
initial amplitude of the fluctuation so as to be as = 0.8 in the ACDM model. 

between the limiting flux and halo mass in the SPT survey [107] is given by 



where cIa(z) is the angular diameter distance. In figure 19, we show the expected redshift 
distributions of clusters in the SPT survey in redshift bins of Az = 0.1 including Poisson 
errors for clusters above the limiting halo mass Mn m , which is determined by eq. (6.32). In the 
left panel, we compare the ACDM model and the kinetic braiding model with n = 5, where 
we fixed O /i 2 = 0.1334, and adopted a$ = 0.8 for the ACDM model and a$ = 0.75, 0.8, and 
0.85, respectively, for the kinetic braiding model. In the right panel, we compare the ACDM 
model and the kinetic braiding model with n = 5 and = 0.27, 0.28, and 0.29, respectively, 
in which we fixed the initial amplitude of the power spectrum for the kinetic braiding model 
so that cr 8 = 0.8 for the ACDM model. 

As demonstrated in figure 19, the redshift distribution of the high-redshift clusters is 
practically useful to distinguish between the ACDM model and the kinetic braiding model, 
depending on the cosmological parameters and the model parameters. However, figure 19 
also demonstrates that the degeneracy between n and ag is problematic. Therefore, in 
order to distinguish between the ACDM model and the kinetic braiding model, we need an 
independent constraint on the amplitude of the perturbations on the kinetic braiding model. 
This will be obtained by combining a measurement of the cosmic microwave background 
anisotropies or a measurement of the cosmic shear statistics from weak lensing surveys. 
The Fisher matrix analysis in section 6.5 estimates a constraint on f^o and n, assuming 
the amplitude of the cosmological perturbation can be determined from those independent 
observations. 




(6.32) 
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Figure 20. The 1-sigma (dashed curve) and 2-sigma (solid curve) contours in the Oo — n plane from 
the Fisher matrix for the cluster distribution assuming like the SPT survey with redshift measurements 
(see section 6.4). The target parameter is n = 5 and iloh 2 = 0.1344. We fixed the spectral index 
n s = 0.96 and the amplitude of the initial cosmological perturbation so that erg = 0.8 in the limit of 
the ACDM. 

6.5 Fisher Matrix Analysis 

In this subsection, we investigate a constraint from the SPT survey predicted with the Fisher 
matrix analysis. The Fisher matrix for the cluster counts (e.g., [109]) can be written as 



where N m is the number of galaxy clusters in the mth redshift bin of the range z m < z < z m +i, 



In the Fisher matrix analysis for the cluster counts, we consider 2 parameters, n and Qq, and 
adopt the target parameter n = 5 and HqH 2 = 0.1344. The other cosmological parameters 
are fixed as those in section 5.3. 

Figure 20 shows the 1-sigma (dashed curve) and 2-sigma (solid curve) confidence con- 
tours in the n and CIq plane, assuming like the SPT survey whose sky coverage is A = 
4 x 10 3 deg 2 and the limiting flux, as explained in the previous section. From figure 20 one 
can see the similar feature as that in figure 16. The constraint on n is not very tight. How- 
ever, the 1-sigma (2-sigma) error in determining n is An ~ 10 (15) when n = 5. Thus, it is 
possible to distinguish between the kinetic braiding model and the ACDM model if n ^ 10. 
However, the constraint on the model parameter n become very weak when the target value 
of n becomes large because the kinetic braiding model approaches the ACDM model. 

7 Conclusion 

In this paper, we studied the cosmological consequences of the kinetic braiding model, whose 
self-interaction term is in the form, G((j), X)0(p oc X n O(j). This model allows the existence of 




(6.33) 




(6.34) 
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a self- accelerating solution without ghost and instability for perturbations. Due to the Vain- 
shtein mechanism, general relativity recovers on small scales, which ensures the consistency 
with the solar system tests. In this model with assuming the attractor condition (2.17), we 
found that the background evolution is identical to that of the Dvali- Turner's model, which 
approaches the ACDM for n equal to infinity. We also found that the linear growth history 
of the cosmological perturbation is identical to that of the ACDM model for n = oo. Using 
the la supernovae data [72] and the WMAP CMB data [4], we obtained a constraint on the 
model parameters of the kinetic braiding model. The kinetic braiding model with n = 1 
requires rather high matter density parameter, which is consistent with the result of [43]. We 
found the model with the higher n can be consistent with the observations. We also found 
that the sound velocity of the galileon field perturbation becomes zero for n = oo. This pre- 
vents the perturbation of the galileon field from propagating, which breaks the quasi-static 
approximation when n is large. 

By solving the full equations of the linear cosmological perturbations numerically, we 
confirmed that the quasi-static approximation can be used for sub-horizon modes k XL 0.01 h 
Mpc -1 as long as n ^ 10. We investigated the evolution of the linear density perturbations 
as well as the spherical collapse in the nonlinear regime of the density perturbations, which 
are important in order to distinguish between the kinetic braiding model and the ACDM 
model when n is small. We also found a useful analytic expression for the growth index, in a 
form expanded in terms of the matter density parameter. The theoretical prediction for the 
large scale structure is confronted with the multipole power spectrum of the luminous red 
galaxy sample of the Sloan Digital Sky survey. The current constraint is not very tight, but 
a future redshift survey like the WFMOS/SuMIRe PFS survey would be potentially useful 
to distinguish between the kinetic braiding model and the ACDM model if n is small. 

Based on the quasi-static approximation and the sub-horizon approximation for model 
n 10, we also investigated the non-linear evolution of density perturbations using the spher- 
ical collapse, where the nonlinear effect of the galileon field perturbation plays an important 
role through the Vainshtein mechanism. We showed that the effective gravitational constant 
in the perturbation including the nonlinear effect approaches the usual Newton's constant af- 
ter the density becomes high. We also investigated a prospect for the number count of galaxy 
clusters in the SPT survey, which will be a useful test by combining other measurements like 
the cosmic microwave background anisotropies or the cosmic shear statistics of weak lensing 
surveys. 

The kinetic braiding model in the present paper reduces to the ACDM model for n = oo, 
which we demonstrated at least as for the background expansion and the linear cosmological 
perturbations. In this meaning, the kinetic gravity braiding can be unbraided by taking a 
large value of n. This is an interesting feature that might be useful to understand the origin 
of the cosmological constant. 
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A Linear Perturbation Equations 



In this appendix, we summarize perturbation equations for the general Lagrangian (2.1). 
The (0, 0) component of the Einstein equations is given by 



Mi 



PI 



-6H($ - H^>) 



K+84 - K X SX - K^dcj) - KxxcfSX 



+ 2Gs5X + 



G x 



+ G^ x r{SX -ZH<t>5<t>) 

12 



30 d $ - + 9H<f 



- SGxxHjfSX - Sp, 



(A.l) 



where 5X is defined by eq.(4.3) and bp is the matter density perturbation. The (0,i) com- 
ponent of the Einstein equation is 



M, 



pi 



2($ - HV) 



(K x - 2G ) 



+ 3H6(/>) +5q, 



(A.2) 



where 5q describes the velocity field of the matter component (See below for details). The 
traceless part of the component of the Einstein equation gives (4.5), while the trace 

part yields 



M 



pi 



1 



2(3H 2 + 2H)V + 2H^ - 2$ - 6H<5> + -^V 2 (^f + $) 



+ K X SX - 2G 4l 5X 



+ G X 



-Gxx<P 2 4>5X + 5p 



(A.3) 
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where bp is the perturbation of pressure. The perturbation equation for the galileon field is 
given by 



(2G - K x ) 



3(^6 - 0* - 2(6 + 3H<j>W + 56 + 3H86 - 4rV 2 c 



+ (K H - K HX $ + G H ^ 2 )56 - (Kxxx + 3G X xxH^)4> 2 '(p5X 
- {K^x - 2G H ) \SX + (4> + W4>)86] - %^ 2 (<5I + ," 



+ K XX 



+ 3H<j))^ 2 ^ - 4> 2 56 - 3(4> + H<j>)<j)56 



G x 



- 12 j(# + 3H 2 )6 2 + 2H66^ - ^ V 2 ^ + 6iT#0 
+ Q^H '6' + (ij + 3# 2 )<4<ty - ^(0 + 2F(/»)V 5 

— Gfix 



3<j) 2 <5> + 6(0 + 3H<p)<p<5> - 9H<j) 2 ty 

12 



(3$ + ¥) + 6(6 - H6)6 2 ^ - 6 2 56 - (46 - 3H6)656 

12 



+ 3< 2H66 + (H + 3H 



2\i2 



Gxx 



30 3 0$ - 3H<j)^ - 3\ SH^'6' + (H + 3H 



r2\I4 



+ 3H6 3 56 + 3<^ 5H6 2 6 + (H + 3H 



2\13 



+ G^xft 



5X + (6- 3H6)56 



+ G rj>X x6 2 (6-3H6)5X -3H6656 =0.(A.4) 



Let us summarize the perturbations of the matter component, whose evolution equation 
is given by the energy-momentum conservation V^T^ = 0. The perturbation of the energy- 
momentum tensor is given by 5T$ = —5p, 5T® = di5q, and 5T l - = 5p5j, respectively. The 
energy-momentum conservation yields 



bp + 3H(5p + 5p) + 3$(p + p) + -^5q = 0, 

or 

6q + 3H5q + (p + p)* + 5p = 0. 



(A.5) 
(A.6) 



For the matter perturbations, we can set 5p = 0. It is convenient to define the quantity 
5q = —pv. In terms of the density contrast 5 = 5p/p and v, these equations can be written 

as 



• V 2 
5 + 3<S> T v = 0, 

ar 

v = 



(A.7) 
(A.8) 



Taking the time derivative of (A.8) and using (A.8), we obtain the evolution equation for the 
matter density perturbation, 



. V 2 

5 + 2HS = - 3$ - 6H<S>. 

a 1 



(A.9) 
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Finally, we introduce the gauge-invariant density perturbation, A c = 6 + 3Hv, and the new 
variable, B = — <E> + Hv. Then, the matter perturbation obeys the following equation, 

V 2 

A c + 2HA C = + W + 6HB. (A.10) 

er 
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